Perturbative treatment of the multichannel interacting resonant level model in steady 

state non-equilibrium 
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We consider the steady state non-equilibrium physics of the muhichannel interacting resonant level 
model in the weak coupling regime. By using the scattering state method we show in agreement 
with the rate equations that the negative diflerential conductance at large enough voltages is due 
to the renormalization of the hopping amplitude thus of the vertex. 
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Since the experimental realization of quantum dots the 
electronic transport in mesoscopic systems has attracted 
fastly growing interest^. In these systems the short range 
Coulomb interaction plays an important role and may 
result in strong correlations between the electrons on the 
dot and the electronic leads. 

In the following we consider a single level quantum dot 
with spinless electrons thus only the interaction between 
the dot and the leads is taken into account. The dot is 
attached to several lead electrodes thus hopping between 
the dot level and the neighboring sites of the leads are 
present in addition to the short range Coulomb interac- 
tion acting between the dot electron and the electrons on 
the neighboring sites of the leads as sketched in FiglT] 




FIG. 1: Possible realization of the multichannel interacting 
resonant level model. A single level quantum dot is attached 
to several lead electrodes. Hopping between the dot level 
and the neighboring sites of the leads A'' = 1 , 2 are present 
in addition to the short range Coulomb interaction acting 
between the dot electron and the electrons on the neighboring 
sites of all A'' electrodes. The voltage is applied between leads 
1 and 2. 



Without losing generality, we can assume that the hop- 
ping is restricted to the leads 1 and 2 while the Coulomb 
interaction acts to all of the N leads. The model thus 



described by the Hamiltonian 

k,i=l...N 



k,i=l,2 



k.i=l...N 



(1) 



where cj^^ i'^ki) creates (annihilates) an electron in the 
lead i with momentum k, d stands for the dot electron 
annihilation operator, to is the hybridization amplitude 
between the dot and the leads 1 and 2, while U is the 
strength of the Coulomb interaction acting between the 
dot and all the leads. The substraction of 1/2 in the 
last term is due to practical reasons: in this form the 
Hamiltonian obeys electron-hole symmetry if the local 
level is tuned to Ed = 0. 

The equilibrium physics of the model has been inten- 
sively studied in the recent decades^i^i^iS . The steady 
state non-equilibrium situation -beyond perturbation 
theor y^'^ has been studied by exact and numerically ex- 
act approaches including Bethe Ansatai, conformal field 
theory and time dependent density matrix renormaliza- 
tion group^iSiiiS. 

By applying these methods the I-V characteristics of 
the model has been calculated where the external voltage 
applied between the conduction electron leads labelled 
by 1 and 2 (left /right). Even if the exact methods give 
correct answers it is hard to get insight into the phys- 
ical mechanisms and the different competing processes. 
The perturbative methods can provide help in the bet- 
ter understanding^'S.. Until now it was applied only to 
the relatively weak voltage bias case. Recently, at high 
bias a rather surprising phenomenon, namely a nega- 
tive differential conductance has been suggested by exact 
methods^iii. The aim of the present paper is to provide 
hints by perturbative methods how that negative differ- 
ential conductance is established. 

Following the earlier work the scattering state method 
combined with time dependent perturbation theory will 
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be applied. That method is suitable for out of equilib- 
rium situations as in the method any initial state can be 
considered. The renormalization is taken into account by 
a generalization of Anderson's poor man's scaling includ- 
ing the corrections in ncxt-to-lcading logarithmic order 
and the renormalization of the initial and final states^^. 

In the spinless model two competing phenomenon oc- 
cur. 
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(i) If the impurity level is occupied then the repul- 
sive Coulomb interaction pushes away the electrons 
in the leads from the vicinity of the dot therefore 
creating more unoccupied states which helps the 
dot electron to tunnel to one of the leads. In this 
way the Coulomb interaction enhances the hopping 
rate. 



FIG. 2: The bare Coulomb vertex (a) and hybridization ver- 
tex (b). The solid lines represent the conduction electron 
propagators in different channels and the dashed line the elec- 
tron on the dot. Diagram (c) is the self energy, (d) the vertex 
correction to the hopping in first order and (e) in second or- 
der. The flow of time agrees with the direction of the dotted 
line. 



(ii) When the occupation of the dot is changed then 
it has to be followed by the rearrangement of the 
electrons in the leads due to the Coulomb interac- 
tion. This rearrangement is similar to Anderson's 
orthogonality catastrophe takes a long time to be 
completed. This mechanism essentially reduces the 
hopping rate. 

These two processes are competing in case of repulsive 
Coulomb interaction U. In perturbation theory the for- 
mer one appears already in first order of U, while the 
latter one in second order only. Therefore a crossover 
is expected by increasing C/i^ In case of two leads that 
crossover occurs in the medium/strong coupling regime 
which is already outside of the validity region of the per- 
turbative treatment. The crossover can be pushed down 
to the weak coupling regime by increasing the effect of the 
Coulomb repulsion by increasing the number of screening 
channels (leads) e.g. to iV = 10. 

As in Ref,5 the self energy and the vertex corrections 
are given by the skeleton diagrams in Fig[2]The solid lines 
represent the conduction electron propagators in different 
channels and the dashed line the electron on the dot. 
Diagrams (a) and (b) are the bare vertices (c) is the self 
energy, (d) the vertex correction to the hopping in first 
order and (e) in second order. The flow of time agrees 
with the direction of the dotted line. 

The self energy and vertex corrections show different 
features concerning the role of the applied voltage. In 
the self energy diagrams the occupation of the dot is not 
changed thus the value of the applied voltage is irrelevant 
as only electron-hole pairs are created at one of the Fermi 
surfaces. In the vertex correction the occupation of the 
dot is changed by a hopping event thus the voltage V is 
relevant and occurs as an infrared cutoff. 

The vertex correction in first order of U (see Figl5](d)) 
is similar to a Hartree diagram thus it is independent of 
the incoming energy to thus it is constant but depends 
on the voltage as Ut\og[D/{eV + uJc)] where uJc is the 
infrared cutoff arising from the width of the resonance-'^. 
It is not necessarily the case for the correction of higher 



order. E.g. the correction depicted in Figl^jje) gives 
Uho log log 

eV \0J\ + UJr 



(2) 



for uj ^ eV. That may lead to difficulties to give the 
renormalization of the vertex for large energies. However, 
in the transport only the electrons in the energy window 
—eV/2 < oj < eV/2 play a role thus the appearance of 
such energy dependence is out of the focus of our interest. 
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FIG. 3: The normalized impurity density of states for = 10 
channels at goU = 0.1 and different values of the bias voltage. 



r = 2-KQotQ. 



The scaling equations provide the expressioi 
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t(w)-to — ^ , (3) 
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where qq is the conduction electron density of states and 
it is assumed that the width of the resonance F can be 
neglected in comparison with the applied voltage eV . 

As in Ref[^ the rate equations are used to determine 
the steady state current between the dot and one of the 
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FIG. A: I — V characteristics with the clear evidence of the 
non-monotonic behavior for the A*' = 10 channel interacting 
resonant level model for different values of the Coulomb re- 
pulsion. 



leads due to the voltage eV applied on the leads (see Eqs. 
(28)-(31) of RefH). In the scattering state method it is 
more feasible to calculate the scattering amplitude of an 
electron going from one of the leads to the other one. In 
the framework of the time ordered diagrams the first and 
last hopping must be picked up and the propagation of 
the electron on the dot can be described by the inverse 
life time V{ljj)/2 which contains the contribution to all 
orders in the hopping to the leads. In the initial state 
the impurity can either be occupied or unoccupied but 
the results are very similar. Thus that amplitude is 
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uj + iT{ijj)/2 



(4) 



and in this way the current is 



/(el/) 



eV/2 
~eV/2 



uj + iT{uj)/2 



duj 
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where the current flows in the energy window —eV/2 < 
u! < eV/2 if the voltage bias applied symmetrically. 
t{uj) appearing in ([5]) is determined self consistently, i.e. 
its RG flow is terminated either by w, the voltage eV 
or by r(a;) itself according to ([3]). Using the relation 
2Trgot'^{uj) — T{ll>) and 
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r(w)/2 



7rw2 + (r(^)/2)2 



eV/2 



I{eV) = 27re 



(6) 
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Comparing Eqs.jl]) and ([6]) it is interesting to notice that 
a factor is now incorporated by the density of states 
Qd- The result given by ([6]) is just the current going to the 
dot from e.g. the left electrode. That holds only in the 
special case with L/R symmetry (t^ = tu). It is worth 
mentioning that in this special case the results obtained 
by using rate equations and the scattering state method 
coincide. 

Considering the voltage dependence of the current two 
factors play role. According to Eq.Q the voltage de- 
pendence of the hopping rate i(a;) is determined by the 
vertex correction which is reduced if the voltage is in- 
creased. On the other hand the renormalized density 
of states Qd is increasing with increasing voltage in the 
relevant energy window (see Fig [3]). As the final result 
shows current reduction at large enough voltages, there- 
fore the dominating process is the vertex renormalization 
(see FiglH). 

In conclusion, we have shown that remarkable physics 
of negative differential conductance in the interacting res- 
onant level model previously studied by exact methods 
can be understood on the ground of perturbation theory. 
Such a treatment of the problem -though it is restricted 
to the weak coupling regime of the model- permits us to 
gain insight into the mechanisms leading to negative dif- 
ferential conductance at large voltages. Our conclusion 
is that the underlying mechanism is the vertex renormal- 
ization, therefore, the mechanism leading to a negative 
differential conductance takes its origin in the equilib- 
rium dynamics of the model. As the interacting resonant 
level model has become a benchmark model for different 
theoretical methods dealing with non-equilibrium physics 
of quantum impurity problems, it would be nice to check 
whether the more precise RG method o^^i^^'^^i^^'^^ pro- 
vide conclusion similar to ours. 
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